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DYNA:  DYNAMIC  STORAGE  ALLOCATION 

IN  FORTRAN  FOR  THE  IBM/ 360  OPERATING  SYSTEM 

ABSTRACT 

This  report  describes  an  I.B.M./360  Assembly  Language  routine  to 
allow  for  dynamic  allocation  of  core  storage  in  Fortran  programs.  The 
report  describes  in  detail  the  implementation  and  use  of  the  dynamic 
allocation  subroutine  DYNA. 

A  fairly  detailed  knowledge  of  FORTRAN  is  assumed  in  the  report. 
The  dynamic  routine  may  be  implemented  by  the  FORTRAN  programmer 
without  a  detailed  underatandlng  of  the  Assembly  Language  routine 
DYNA.  A  complete  description  of  DYNA  is  included.  The  discussion 
is  technical  and  assumes  a  knowledge  of  I.B.M./360  Operating  System 
'lacros  and  Assembly  Language. 
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DYNA:  DYNAMIC  STORAGE  ALLOCATION 
IN  FORTRAN  FOR  THE  IBM/ 360  OPERATING  SYSTEM 

Alan  C.  H.  Kao 
and 

Charles  F.  Wall 
University  of  Hawaii 

1.  INTRODUCTION 

A  major  limitation  of  Fortran  is  that  storage  for  data  is 
allocated  at  compile- time  by  the  programmer.  Essentially,  some 
oplmum  dimensions  are  selected,  the  program  is  compiled  and  fixed 
amount  of  storage  is  made  available.  If  a  particular  problem  requires 
less  space  than  is  allocated,  storage  is  wasted.  If  more  space  is 
required  than  has  been  alloc -.ted,  the  program  must  be  redimensioned 
and  recompiled.  In  today's  multlprogrammed  systems,  either  situation 
is  an  inefficient  use  of  system  resources. 

A  solution  to  the  problem  is  to  allocate  storage  for  data  at 
execution  time,  dynamically.  The  general  rationale  is  that  allocation 
of  storage  for  data  is  a  dynamic  function  that  varies  with  each  use  of 
a  particular  program. 

2.  DYNA:  Dynamic  Storage  Allocation  in  Fortran  for  the  IBM/360 

DYNA  is  an  assembly  language  subroutine  that  provides  an  inter¬ 
face  between  a  Fortran  program  or  subroutine  and  the  operating  system' 
(O.S.)  main  storage  manager.  Each  time  DYNA  is  called,  storage  is 
allocated  and  made  available  to  the  Fortran  subroutine. 


2. 


The  next  four  sections  discuss  ;he  O.S.  Getmain/Freemain  facil¬ 
ity,  the  interface  for  these  functions,  and  the  Error  Conditions  en¬ 
countered.^ 


2.1  Fortran  Convention 

The  standard  Fortran  call  ia  used  to  call  DYNA  and  to  allocate 
storage.  The  order  of  the  parameters  is  as  follows: 

CALL  DYNA  (SUB,  ARRAY1,  NSIZE1,  ...»  ARRAYN,  NSIZEN) 


where 

SUB  is  the  name  of  the  Fortran  subroutine  the  storage  is 
allocated  to 

ARRAYN  la  the  Nth  array  to  be  allocated 

NSIZEN  ia  the  count  of  storage  in  words  of  ARRAYN. 

DYNA  is  called  by  the  calling  Fortran  routine  which  in  turn  calls 

the  Fortran  subroutine  SUB.  SUB  must  be  declared  EXTERNAL  in  the  calling 

routine  and  the  parameters  must  not  appear  in  COMMON. 

The  parameter  list  is  a  list  of  the  addresses  of  epch-parameter 
in  the  subroutine  call.  The  first  parameter  is  the  address  of  the  sub¬ 
routine  to  be  called  by  DYNA.  The  last  parameter  has  its  high  order 
byte  equal  to  128.  The  call  to  DYNA  consists  of  loading  the  address 
of  the  parameter  list  in  general  register  1,  loading  the  address  of 
DYNA  in  register  IS  and  then  branching  to  register  IS  and  linking 
register  14.  Register  13  points  to  a  save  area  within  the  calling 
routine  (standard  O.S.  linkage  convention). 

See  References. 
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2 . 2  0  .S_,  Getmain/Freemain 

Getmain  la  an  SVC  which  requests  the  control  program  to  assign 
ownership  of  an  area  of  main  storage  to  a  routine.  This  SVC  has  many 
options  and  Is  conveniently  coded  as  a  macro  in  assembly  language. 
However,  It  has  been  hand  coded  in  DYNA  using  the  R  and  VU  forma.  The 
R-form  requests  a  specific  amount  of  main  storage  in  register  0.  The 
SVC  returns  the  address  of  this  area  in  general  register  1.  The  VU- 
form  requests  an  amount  of  main  storage  between  two  ranges.  If  the 
request  cannot  be  satisfied,  the  control  program  will  terminate  the 
calling  routine  (see  Error  Conditions).  For  this  form,  register  1 
points  to  a  parameter  list  with  the  appropriate  flag  specifying  the 
VU-request . 

The  Freemaln  SVC  frees  the  storage  assigned  by  the  Getmain. 

Both  the  R  and  V  forms  are  used  and  are  hand  coded. 

2.3  The  Interface 

Upon  entry  to  DYNA,  the  general  registers  are  saved  in  an  area 
provided  by  the  calling  routine.  Then  the  Getmain  SVC  issued  to  obtain 
a  work  area  is  appropriately  Initialized. 

The  next  phase  is  to  process  the  parameter  list  and  save  the 
address  of  the  calling  routine.  The  parameter  requests  are  rounded  up 
in  preparation  for  the  double-word  alignment.  Once  the  total  size  of  the 
arrays  hus  been  calculated,  the  second  Getmain  is  issued.  The  main 
storage  block  allocated  is  used  by  DYNA  to  reallocate  array  addresses. 

The  original  array  addresses  are  replaced  with  the  new  double-word 
aligned  addresses.  Thus,  the  initial  parameters  corresponding  to  the 


4. 

addresses  passed  to  DYNA  of  the  arrays  are  replaced  by  the  addresses 
of  the  main  storage  area  to  be  used. 

DYNA  then  calls  the  subroutine  whose  address  is  given  by  the  first 
parameter  In  the  list  passed  to  DYNA.  On  return  from  the  called  sub¬ 
routine  DYNA  frees  the  main  storage  allocated  for  the  arrays  and  re¬ 
establishes  the  linkage  to  the  calling  program.  DYNA  then  frees  the 
work  area  allocated  and  returns  to  the  calling  routine. 

Since  DYNA  obtains  a  work  area  (activation  list)  upon  entry  and 

frees  this  work  at  exit,  DYNA  Is  a  re-entrant  routine  and  may  be  called 

2 

by  more  than  one  program  In  any  order. 

2.4  Error  Conditions 

If  there  Is  no  main  storage  available  the  system  will  abnormally 
end  (ABEND)  with  a  code  of  S80A  or  S804.  Other  codes  are  possible  but 
have  yet  to  be  encountered. 

If  there  Is  some  storage  available  but  not  enough  to  satisfy 
the  request,  DYNA  calls  $ERR1,  an  error  routine  written  In  Fortran 
which  has  as  an  argument  a  vector  with  the  following  values:  the 
address  of  the  core  that  was  allocated,  the  amount  the  system  could 
allocate,  the  minimum  requested  size  (usually  zero)  and  the  maximum  size 
requested. 

$EUR2  Is  called  if  the  parameter  list  has  the  Improper  number 
of  arguments.  There  should  always  be  an  odd  number  of  arguments  since 
the  first  argument  is  the  name  of  the  Fortran  subroutine  called  by 
DYNA  and  all  other  arguments  are  the  paired  names  and  sizes.  This 

2 

See  References. 
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error  exit  occurs  In  the  debugging  stage  and  does  not  return  control. 

A  message  is  printed  and  control  is  returned  tc  the  system. 

3.  PROGRAMMING  DYNAMIC  ROUTINES 

The  use  of  DYNA  is  relatively  straightforward.  The  subroutine 
(DYNA)  is  reentrant  and  may  be  nested  to  any  level.  The  tradeoff  is  in 
the  time  required  to  call  DYNA  in  comparison  to  its  efficient  use  of 
space. 


3.1  Some  Rules 

As  mentioned  in  2.1  standard  Fortran  is  used  to  call  DYNA. 

The  name  of  the  subroutine  to  which  DYNA  will  pass  control  must  be 
declared  EXTERNAL  in  the  routine  calling  DYNA.  The  narameters  in  the 
call  statement  must  not  appear  in  COMMON.  If  an  array  is  to  be  DOUBLE 
PRECISION,  the  size  parameter  must  reflect  this  (e.g.  ten  words  DOUBLE 
PRECISION  must  be  passed  as  twenty  words) .  DYNA  always  allocates  on  double- 
word  boundaries.  If  a  particular  array  is  doubly  subscripted  the 
following  form  is  not  allowed  in  the  called  subroutine: 

DIMENSION  X(l,l) 

Instead,  the  appropriate  dimensions  should  be  passed  in  COMMON  and  used 
in  the  called  program  in  the  following  manner: 

DIMENSION  X(N,M) 

COMMON  N,M 


For  vectors,  either  form  is  acceptable. 
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3.2  An  Added  Advantage 

If  your  program  follows  the  normal  Fortran  convention  of  storing 
arrays  columnwise,  the  arrays  may  be  dcubly-subscrlpted  and  calls 
to  subroutines,  as  In  the  Scientific  Subroutine  Package,  may  be 
called  without  conversion.  Since  we  always  allocate  the  exact  amount 
of  storage  needed  by  each  array,  the  array  may  also  be  treated  (In 
separate  routines)  as  a  vector.  Since  it  Is  useful  to  use  doubly-* 
subscripted  arrays  for  input/output  and  certain  computations  and 
singly-subscripted  arrays  (vectors)  for  other  computations  It  Is 
useful  to  be  able  to  eliminate  conversion  problems.  (For  a  more 
complete  discussion  see  IBM/Scientif i~.  Subroutine  Package-Programmers 
Guide) . 


3.3  A  Sample  Program 

Consider  a  program  which  calculate  means,  standard  deviations  and 
correlations.  A  simple  call  to  the  IBM/SC?  subroutine  C0RRE  will  suffice. 
All  1/0  will  not  be  shown  as  the  purpose  is  to  demonstrate  the  use  of 
DYNA.  We  shall  consider  that  we  are  calling  the  DOUBLE  PRECISION  version 
of  CORRE. 

C  —  THIS  IS  THE  MAIN  ROUTINE  — 

COMMON  N,  M 
EXTERNAL  SUBC 

C  —  READ  C0NTR0L  CARD  AND  ALLOCATE  C0RE 
READ  (5,  10)  NR0WX,  NC0LX 
10  FORMAT  (213) 

N  -  NROWX 


M  -  NCOLX 
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N1  -  NROWX  *  2 
N2  -  NC0LX  *  2 
N3  -  N1  *  N2 
N4  -  (N2  *  (N2  +  l))/2 
C  —  CALL  SUBC  VIA  DYNA 

CALL  DYNA(SUBC,X,N3,XBAR,N2,STD,N2,RX,N3,R,N4,B,N2,D,N2,T,N2) 

ST0P 

END 

C  —  THIS  IS  SUBROUTINE  SUBC 

SUBROUTINE  SUBC(X,NA,XBAR,NB,STD,NC,RX,ND,R,NE,B,NF,D,NG,T,NH) 
C  —  NOTICE  THa.*  IN  THE  SUBROUTINE  WE 
C  —  CANNOT  HAVE  DUPLICATE  DUMMY  NAMES 

DIMENSI0N  X(N,M) ,XBAR(1) ,STD(1) ,RX(1) ,R(1) ,B(1) ,D(I) ,T(1) 
C0MM0N  H,M 
C  —  READ  IN  DATA 
D0  20  I  -  1,  N 

20  READ  (5,10)  (X(I,J),  J  -  1,  M) 

10  FORMAT  (10F8.0) 

C  —  CALL  CORRE 

CALL  CORRE  <N,M,I,X,XBAR,STD,RX,R,B,D,T) 

C  —  PRINT  RESULTS 


RETURN 

END 


3 . 4  User  Documentation 


Primarily,  the  user  must  know: 

a)  The  formula  to  calculate  the  dynamic  area;  and 

b)  the  amount  of  static  storage  used  by  the  program  and 
buffers. 

A  simple  test  run  with  size  -  1  for  all  dynamic  data  arrays  will  give 
the  base  (e.g.  amount  of  static  storage).  To  be  safe  this  should 
be  rounded  up  to  the  next  even  K  bytes.  The  formula  is  determined  by 
the  particular  program.  For  example,  the  formula  for  our  sample 
program  would  be  given  as  follows: 

SIZE  (in  bytes)  -  (2  *  N3  +  5  *  N2  +  N4)  *  4 

where 

Ml  -  Number  of  rows  in  X  times  2 
N2  -  Number  of  columns  in  X  times  2 
N3  -  N1  *  N2 
N4  -  (N2  *  (N2  +  l))/2 

the  region  is  •  SASE  +  SIZE  rounded  up  to  next  even  K  bytes. 
This  is  the  estimate  the  user  should  supply  on  the  EXEC  card. 

C  —  TO  SUPPLY  LESS. 

C  —  RESULTS  IN  AN  30A; 

C  —  TO  SUPPLY  MORE: 


C  — 


IS  A  WASTE  FOR  SURE: 
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********  ******  **  *  ******  *  ********  *********  *********  *  ******** ****** *** * 

OYNAMIC  CORE  ALLOCATION  PROGRAM  FOR  FORTRAN  PROGRAMS 
_ US_IN£_  OS. .  MACRO  GETMAIN  FOR  THE  ^ALLOCATION  MECHANISM, _ 

2.  DOUBLE-WORD  BOUNDARIES 

CO'RE*  GUESTIMAT  ION  ALGOR  I THM  * 


LINKAGE  EOITOR  MAP_  TOTAL  LENGTH _ 

SYSIN/SYSPRINT  4K  BYTES 

FOR  N  OD  CARDS  N*2K  BYTES 

SUM  ( BLKS I ZE*BUFNQ  )  M 

TOTAL  REQUIREMENTS  (TOTAL  LENGTH) +4K+N*2K+M 


_***«*j|t****  ************************************************  ************ 

DYNA  C  SECT 

STM  14,12*12(13)  SAVE  THE  GENERAL  REGISTERS 

_ LR  12,15  _  _USE  GR  12  AS  THE  BASE  REGISTER 

USING  OYNA, 12  TELL  THE  ASSEMBLER  THE  BASE  IS 

LR  2,1  SAVE  REGISTER  1  FOR  LATER. 

_ LA  _0,  SAVE  SI  ZE _ GR  0  HAS  THE  SIZE  OF  SAVE  AREA 

8AL  1 ,*+4  *  IBM  CONVENTION  FOR  GETMAIN  (R) 

SVC  10  GET  MAIN  SVC 

_  USING  _  OYNASAVE.il  _ 

~LR  11,1  . USE  GR'  11  AS  SAVE  BASE  REGISTER 

LR  1,2  RESTORE  GR  1  TO  ORIGINAL  VALUE 

_ ST  13.SAVEAREA+4  _  SAVE  GR  13  FOR  LINKAGE 

LA'  •  "  13, SAVEARE A  RESET  THE  LINKAGE 

ST  1 , SAVEGR1 

*  * 

*  INITIALIZATION  OF  THE  PARAMETER  LISTS,  AND  CONSTANTS  . . 

* 

LA  _ 3  .MINSIZE 

ST  ""  3  ,GET  SIZE  "  " 

LA  3 , FREEAREA 

_ ST  3.GETAREA 

ST  3,  ERR  LIST  "  "  “ 

MV I  GETMODE, X*  EO* 

MV I  GETSP,X»00» 

"XC  FREEAREA(16)  .FREEAREA .  **"' 

EJECT 

********************************* 4* ***************************** ****** 


*  NOW  START  ALLOCATING  THE  CORRECT  ADDRESSESS. 

* 

********************************************************************** 

L  5.SAVEGR1 

L  4,0(0, 5)  GET  ADDRESS  OF  REAL  MAIM  PROGRAM 

'  L  4, 0(0, 4)  GR  4  HAS  THE  ADDRESS  OF  REAL  ONE 

ST  4.SAVEPAIN 

LA  ‘  5,4(0, 5)  GR  5  POINTS  TO  THE  ARRAY  ADDRESS 

$T  -  -  --5fSAVEGRl  *  -  SAVE  THE  NEW  PARML1  ST  PNTR. 

* 

*  WANDER  THRU  THE  PARAMETER  LIST  TO  DETERMINE  THE  MAX  CORE  REQ'D 

*''  '  ■“  . . . 

LA  0,0  ZERO  THE  WORK  REGISTER. 

LA  4,0  ZERO  THE  WORK  REGISTER. 

'  'LOCATES  EQU  . *  "  '  *  -  - - 
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CL! 

0(5),X»80» 

CHECK  FOR  ERRONEOUS  PARM  LIST 

6~E 

L 

ERROR  1 

4, 4(0, 5) 

GET  THE  NUMBER  OF  WORDS  REQ'O 

L 

4, 0(0, 4) 

THE  ACTUAL  NUMBER 

S*RL 

4.1 

NUMBER  OF  WORDS/2 

LA 

4, 1(0, 4) 

SAFETY  FACTOR 

SLL 

4,3 

PREPARE  FOR  DOUBLE  WORDS 

"AR  0,4 . Sum  up  the  number  of  double  worc 

CL!  4(b), X  *  80  *  CHECK  FOR  END  OF  LIST 

LA  5, 8(0, 5)  STEP  FOR  THE  EXT  SET 


LOCATES 

w 

* 

FORCING 

DOUBLE  WORO  ALIGNMENT  FOR  DOUBLE  PRECISION  ARRAYS. 

w 

ST 

0 , MAX  SI ZE 

FOR  THE  GET  MAIN  ROUTINE 

LA 

1 «GET  LIST 

“SVC 

4 . .  ' 

GETMAIN  SVC  CODE 

L 

2, FREEAREA 

Gk  2  HAS  THE  ADDRESS  OF  FREECORE 

L 

3,FREESIZE 

GR  3  HAS  THE  AMOUNT  FO  FREE  CORE 

C 

3 • MAXS  IZ E 

TEST  IS  ALLOCATED  EMOUGH 

BL 

ERROR 

L 

5.SAVEGR1 

RESET  GR  1 

LOOP 

EQU 

* 

ST 

2,0(0, 5) 

CHANGE  THE  OLD  ARRAY  ADDRESS 

L 

6, 4(0,5) 

PICK  UP  THE  ADDRESS  OF  «BYTES 

L  . . 

"6, 0(0, 6)  *'*■ 

GR6  «  NUMBER  OF  WORDS  IN  ARRAY 

SRL 

6,1 

■WORDS/2 

LA 

6, 1(0, 6) 

*W0RD»=W0RDS+1 

SLL*  ' 

6,3- 

*»WORDS=MULTIPLE  OF  8 

AR 

2.6  • 

GR  2  ->  NEW  FREE  AREA 

CLI 

4(5),X»80» 

CHECK  IF  LAST  ARGUMENT 

*"*  *  **“  ~  *** 

“LA  " 

5  *8  (0*  5)  ” 

.  GR  5  ->  NEW  ARRAY  LIN  LIST 

3NE 

LOOP 

EJECT 

* 

* 

NOW  CALL 

THE  REAL  MAIN 

PROGRAM  WITH  THE  NEW  PARAMETER  LIST 

♦ 

. L  . 1 ,  SAVEGR1  - 

L  15.SAVEMAIN 

BALR  1*4, 15  GO  TO  IT. 

*  '  *  . . .  . 

♦  FIRST  FREE  THE  ALLOCATED  CORE. 

*  _  RETURN  TO  THE  MAIN  CALLER  NORMALLY. 

# 

RETRUNS  EQU  * 

L  0,FREE  SIZE 

*  L  .  1 , FREEAREA 

SVC  10  FREE  THE  ALLOCATED  CORE. 

RETURNQ  EQU  _  * 

. -  -l3;SAVEAREA4.4  •  —  -  •  • 

LA  •  0 , SAVESI ZE  FREE  UP  THE  SAVE  AREA 

LR  1,11  RETURN  THE  ALLOCATED  ADDRESS 

$VC  .  10  FREE  MAIN  . ' 

LM  14,12,12(13) 

LA  15,0 

_Br  14 

* 

*  _  _NOT  ENOUGH  CORE  IS  ALLOCATED  FOR  THIS  JOB 

*' . 
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ERROR  EQU  * 

LA  “  ltERRLIST  " 

L  .  1 5»«V ( SERR 1 ) 

_ LA  _  1A.RETRUNS 

BR  15'  " 

¥ 

*  _ INCORRECT  SETUP  OF_  THE  PARAMETER  LIST,  MISSING  »OF  BYTES 

• 

ERRORI  EQU  * 

_ LA _  _14,RETURNQ 

L  15  »*V ( SERR2) 

BR  15 

_  EJECT 

***************  ******************************************************* 


DYNAMIC  SAVE  AREA  AND  CONSTANTS  FOR  PROGRAM 

THE  ADDRESSESES  AND  CONSTANTS  ARE  FILLED  IN  AT  EXECUTION  TIME 


*********  ***********************  ************************************* 


LTORG 

_ 

. . .  '  " ' 

DYNASAVE 

DSECT 

•  f 

t 

SAVEBEG 

EQU 

*  ■ 

DS 

OF 

GETLIST 

EQU 

* 

GETSIZE 

DS 

ALA (M I NS  I ZE) 

PARAMETER  LIST  FOR  GETMAIN 

GETAREA 

DS 

ALA (FREEAREA) 

GETMODE 

DS 

XLl'EO' 

VC  MODE 

GETSP 

OS 

XLl'OO* 

SP  0 

* 

DS 

OF 

FREEAREA 

DS 

ALA(O) 

ADDRESS 

OF  FREEE  CORE  AREA 

FREESIZE 

DS 

1 F '  0* 

NUMBER  i 

OF  BYTES  ALLOCATED 

MINSI ZE 

DS 

IF  *  O' 

MINIMUM 

CORE  REQUIREMENT 

MAXSIZE 

£ 

DS _ 

lF'O' 

MAXIMUM 

CORE  REQUIREMENT 

ERRLIST 

DS 

ALA ( FREEAREA ) 

SAVEAREA 

DS 

1 8F • 0 ' 

SAVEGR1 

DS 

0 IF  *0 • 

SAVEMAIN 

DS 

l  F '  0* 

SAVE END 

EQU 

*  ‘ 

- 

SAVESIZE 

EQU 

SAVEEND-SAVEBEG 

END 
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SUBROUTINE  $ERR 1 ( I  CORE )  SERR1010 

DIME  NS  ION  I  CORE  (4)  '  ~  *ERR1020 

PRINT  10»ICORE( 1 ) • ICORE ( 2) » I  COR  E( 4)  JERR1030 

FORMAT ( ' OOYNA  ERROR:  INSUFFICIENT  CORE  FOR  PROGRAM  '//  $£RR1040 

1  *  "ALLOCATED  CORE  AT  ' «  Z8  »•  OF  •  ,  16,*  BYTES.  BUT  RE  OU I  RED"  *  JERRI050 

2  .16,'  BYTES. '//I  SERR1060 

RETURN  •'  JERR10  70 

END  ~  • '  Tr"'  "  "  ‘  ~  .  SERR1080 


SERR20 10 
SERR2020 
SERR2030 
SERR2040 
JERR2050 


SUBROUTINE  SERR2  . . 

PRINT  10 

FORMAT! 'ODYNA  ERRORJ  INCCRRECT_PARAMET_ER_  LJ^ST  SETUP') 

STOP 

END 
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